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Abstract 

The classical characterization of wave propagation , as a typical con- 
cept for far field phenomena, has been successfully applied to many wave 
phenomena in past decades. The recent reports of superluminal tunnelling 
times and negative group velocities challenged this concept. A new lo- 
cal approach for the definition of wave velocities avoiding these difficulties 
while including the classical definitions as particular cases is proposed here. 
This generalisation of the conventional non-local approach can be applied 
to arbitrary wave forms and propagation media. Some applications of the 
formalism are presented and basic properties of the concept are summa- 
rized. 



Introduction 



The velocities conventionally used for the characterization of wave propagation 
are the phase velocity v p h := u/k, the group velocity v gr := du/dk, the front 
velocity Vf := lim^^Vph, the signal velocity v s with Vf > v s > v gr , the energy 
velocity ve '■= P/p and the phase time velocity v v := x/ (d<p/du), where u denotes 
the angular frequency, k - the wave number, 0- a phase shift accumulated in course 
of propagation, x- the distance of propagation, P - the energy current density 
and p - the energy density. 

*e-mail: drosdow@uni-koblenz.dc 
^e-mail: alfons@uni-koblcnz.de 
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This classical concept [1, 2] and its improvements [3, 4] have been successfully 
applied to many wave phenomena in past decades. The recent reports of super- 
luminal and even negative velocities [5, 6] the advent of near field optics and the 
generation of ultra-short pulses [7, 8], for instance, revealed intrinsic difficulties 
and limitations of the classical concept. It is the purpose of the present note to 
resolve these difficulties by a new local approach for the characterization of wave 
velocities containing the classical Ansatz as special case and being applicable to 
arbitrary wave forms and media of propagation. 

The material is organized as follows: the next section contains a formal discus- 
sion of the classical approach pinpointing its intrinsic difficulties and summarizing 
the goals of the new approach. For this purpose it is completely sufficient to focus 
on the notions of phase and group velocity . This prepares the ground for the 
new approach outlined in the section 2 and an elucidation of its properties by 
application to known examples in section 3. 

The behavior of the wave velocities defined here under relativistic transforma- 
tions is outlined in the section 4. In order to pave the way for applications of the 
formalism to experimental results like those mentioned above, the discussion is 
concluded by a practical interpretation of the definitions and the extension from 
a local to a global evaluation of signal velocities . 

2 Some remarks on wave velocities 

The classical concept of wave velocities analyses wave propagation in the far 
field, where the wave is far away from the source. The wave is described via 
propagating harmonic functions containing a periodic factor of the form e t<p . The 
argument of this periodic functions is interpreted as a phase and identified with 
(t — x/v) via the free wave equation , where v is the phase velocity [I]. 

A second fundamental concept of characterizing wave propagation is that of 
the "group velocity". This concept is - mathematically speaking - a comprehen- 
sive definition for a specified (linear) superposition of solutions of the free wave 
equation with the same periodicity properties usually expressed by the frequency 
distribution of the constituent periodic waves. This definition is also an inher- 
ent far field concept considering "source-free" waves propagating in a strongly 
homogeneous and isotropic medium. This medium is characterized only by its 
"dispersion" (supposed to satisfy additionally the Kramers- Kronig relations), i.e. 
by a dependence of the complex index of refraction n(iv) on the periodicity pa- 
rameter of the wave - the frequency u>. 

The different notions of velocities (and the accompanying notion of a signal) 
developed in this context were a major step forward allowing to accomodate many 
experimental data and to characterize many wave phenomena, (cf. e.g. [2]) It 
is, however, not obvious to what extent these far field concepts can be applied to 
near field problems or any of the other challenges mentioned above. A possible 
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answer to this question simultaneously preserving the classical concepts whenever 
applicable is the focus of the present work. The basic idea is to employ a local 
definition of a propagation velocity . 

Such local approaches can be found in many areas physics and its applications 
in concerned natural sciences. A need for an alternative local approach to the 
analysis of wave propagation was justified already since years in different fields of 
application such as geophysics [21], physical chemistry [20, 25], biophysics [22], 
acoustics [23], dealing with a relative "slow" waves, compared to the optics and 
electronics, thus being long ago in the near-field region, where the actual problems 
require a fine resolution of wave behavior. 

There is thus a justified reason for a local approach to wave velocities since 
this allows to resolve the intrinsic problems accompanying the classical global 
definitions. To elucidate this remark we now list some of these problems. 

The conventional analysis of wave propagation is based on a representation via 
periodic functions by means of Fourier analysis. A well defined phase velocity is 
assigned to each Fourier component equipped with its frequency; a corresponding 
group velocity is subsequently defined for the complete Fourier superposition, i.e. 
for a wave packet. The classical wave theory is based on such definitions relying 
canonically of the phase of periodic functions and on group of such functions . 
This Ansatz is supposed to be the inherent kernel of any characterization of wave 
phenomena and has also been used as a guideline for characterizing general - 
not necessary periodic - electromagnetic pulses [1]. A close inspection, however, 
reveals several shortcomings of this Ansatz with respect to mathematics and 
physics. 

Let us start with the local features of wave propagation . A propagating 
wave is per definition a local space-time distribution which is itself a solution 
of a local differential equation . A definition of a wave velocity as a space-time 
relation should respect this frame being local as well. This is not the case in 
the conventional definitions of wave velocities : these involve the frequency and 
the wave vector, which are (temporally and spatially) nonlocal parameters in 
following the independent space- and time-periodicity. Since these parameters 
do not enter the wave equation , they are usually plugged in from the outside. 
The source free wave equation for instance, admits solutions of arbitrary shapes 
without an a priori assumption of periodicity. 

Simultaneously, the classical approach always takes recourse to a representa- 
tion of an arbitrary solution by a set of periodic harmonic functions , where each 
element of the set does in general not obey the wave equation . Any attempt to 
apply this strategy to waves resulting from non-linear equations such as solitons, 
for instance, elucidates this problem in a clear way since in this case even a linear 
superposition of solutions is not anymore a solution. 

Second, the subject of transmission and velocity of a signal is based on suf- 
ficiently local procedures of measurement [7]. Roughly expressed, the intervals 
between several space-time points are measured. Definitions of velocities based 
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on periodicity parameters can possess generically neither time- nor space-locality. 
Moreover, any Fourier transformation is basically a global object, and all manip- 
ulations concerned are in general mathematically exact only with integration over 
the whole space and whole time, as well as over an infinite frequency band. 

The conventional definitions of the phase and group velocity mentioned above 
[1] are based on the special case of a propagation of periodic harmonic waves in 
homogeneous media. Thus there is already an essential contradiction between 
the local character of wave propagation and the inherent globality contained in 
the definitions of wave velocities. 

Any attempt to construct a local measurable object using Fourier sums or 
integrals contain an essential contradiction as outlined above and provides indeed 
no real locality. This is the source of several problems arising when replacing 
originally local features by " microglobal" ones [3, 12, 13]. 

Even if the approach is supported by mathematical consistency (like the as- 
sumption of an infinite frequency band [1]), it does not lend itself easily to a 
transparent physical interpretation. A nice example thereof is the phenomenon 
of signal transmission in (optically active) media, which triggered a lot of discus- 
sions [3, 11, 10] 

The classical papers [1] and later improvements thereof [13, 4] still contain 
an essential mismatch between local and global wave features, based on several 
mathematically correct but physically misleading definitions (like the notion of 
front and signal velocities [9, 26]). This Ansatz is bound to lead to controversial 
results in applications. For example, it seems not more to be surprising, that the 
subject of "group velocity" fails to describe a propagation of ultra-short pulses [8]. 

These shortcomings of the classical concept of wave velocities motivated the 
Ansatz presented here for the following reasons: 

- The classical harmonic approach represents local wave features in terms of 
generically non-local attributes [15, 16]. Hence: 

- The mathematical equipment is not exact for finite physical values (like a 
limited frequency band); 

- It is therefore not obvious how to apply it to near-field effects, if the size of 
space-time regions are of an order of magnitude smaller than the wave periodicity 
parameter; 

- There is no canonical application thereof for inhomogeneous and anisotropic 
media, as well as for ultra-short pulses of arbitrary form and for nonlinear waves 
[12, 17]; 

- As a consequence, the applicability of this approach for several fields of 
modern quantum optics, nano-optics and photonics is difficult, since these topics 
deal with parameter areas outside its range of validity [5]. Further motivations 
outside of optics have been mentioned above. 

The classical approach is also to some extent problematic from a pure technical 
point of view: it is an essential restriction of generality and universality of the 
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theory enforcing via the Fourier analysis - a number of problematical artifacts. A 
prominent example thereof being the assumption of infinite frequency bands for 
signals. These additional problems are then resolved via the analysis of dispersion 
relations. 

A first step towards a more general concept should be based on an indepen- 
dent alternative approach without doubting the canonical harmonic criteria, i.e. 
leading to the same (verified) results in known cases. 

The aim of the present paper is to establish a transparent criterion for evalu- 
ation of the propagation velocity for a quite arbitrary signal, that does not need 
an explicit representation in terms of periodic or exponential functions at all, and 
with minimal loss of generality in other respects. First of all, we have to recall, 
what is being measured in experiments and what is meant when speaking about 
a "wave velocity", thereby providing the spectrum of wave velocities. 



The following discussion is restricted to a 1 + 1 dimensional space-time for sim- 
plicity. To evaluate the speed of a moving matter point one has to check the 
change of the space coordinate Ax during the time interval At. How can this 
Ansatz be applied to waves ? 

An arbitrary wave is a function defined on the 2-dimensional space-time (x, t) 
and one has no a priori defined fixed points to follow up as in the former case. 
Let us therefore analyze a continuous smooth function ip(x, t) of arbitrary shape. 

In a first step towards a local definition of velocity we consider a vicinity 
U(x) of a certain point x at the certain time t. Further we assume the local 
information about the function ip(x,t) to be measurable, i.e. we should be able 
(at least in principle) to evaluate the values of the function ip(x,t) itself and all 
its derivatives in the point x as well, and in the vicinity U(x). 

A description of propagation is based on monitoring the fate of a certain 
attribute ("labelled point") of this shape, i.e. one finds the same attribute at the 
next moment t\ in the next point X\. 

Let this attribute be labelled by some one-point fixed value ip of the func- 
tion ip(x,t). Suppose, one follows up this local "attribute" and manages a local 
observation of the condition: 



which fixes the space-time points {x, t} where this condition holds (Fig.l). 
Thus we have a function x(t) given in an implicit form in eq.(3.1). 

The first order implicit derivation provides the velocity of this one-point local 



3 N-th order phase velocity (PV) 



ip(x,t) -^o = 
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attribute: 



dx 



«(0)(M):==r = -f|, (3-2) 

dx 



called from here on the "zero-order" or "one-point" phase velocity (0-PV). 

(Here the subject "phase" is treated in sence of local charactreristics in the 
point of shape including all N-order derivatives, N = 0, 1,2, ... i.e. all terms of 
the Taylor expansion in this point). 

It describes in the simplest case the speed of translation of some arbitrary 
pulse, that can be treated as the signal propagation velocity, provided that the 
measured value ip of the amplitude is the signal considered. 

To proceed with the local description of the shape propagation we now con- 
sider as the measurable attribute of tp(x,t) at some time t not the single value 
■00 — V'^oi ^o) a t the point x , but a set of values ip in a local neighborhood 
(like U{x)) close to this point. Then another local attribute can be constructed 
to trace their propagation. If, for instance, one looks at a certain value of the 
first derivative dip(x;to)/dx := xo °f the shape ip, (like dtp(x;to)/dx = 0, i.e. 
at a maximum or minimum point), one should consider the propagation of the 
condition: 

d^(x,t) 

~^x~ = Xo ' (3 ' 3) 
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this provides its propagation velocity via 

d 2 4> 

W(J) (x,f):=-||, (3.4) 

called in view of (3.2) the "first order" or "two-point" phase velocity (1-PV) 
respectively. 

When tracking the propagation of a maximum (minimum), the conditions 
^£ = and < have to hold simultaneously. 

This Ansatz is easily iterated to phase velocities of order N interpreted as the 
propagation velocities of higher order local shape attributes via 

v {N) (x,t) := (3-5) 

leading to N-th order ("N+l-point") PV, describing the propagation of higher 
order local shape attributes. 

The phase velocities so defined are obviously local features depending on 
space-time coordinates. It has to be noted, that any given problem at hand 
might require a particular choice of a PV allowed for by the definitions given 
above. The following examples elucidate these requirements and demonstrate 
the flexibility of these definitions. 

Before proceeding to it should be recalled, that the PV-spectrum has been 
obtained wanting to be able to describe the propagation of an arbitrary pulse in 
terms of local attributes in a medium whose properties depend on several vari- 
ables, especially time and space coordinates being the most prominent and natu- 
ral examples thereof. Any initial shape ip(x ,t ) thus should be deformed during 
the propagation (or evolution, as typicaly encountered for dispersive media). The 
ordinary phase velocity vq therefore is not a relevant criterion to characterize the 
shape propagation and one has to choose an appropriate PV i>(jv)- 

For example, let the pulse ip(x,t ) be subject to damping (Fig. 2). As a con- 
cequence, the zeroth order PV measured in the point x\ gives a magnitude much 
smaller as the same magnitude measured in the point x-i. 

For an amplified signal (Fig. 3, like a signal propagating in a laser excited 
medium), by comparison, the zero order PV from the point X\ provides alto- 
gether even a backward propagation. In both cases an appropriate approach 
would be to apply the first order PV i>(/) which describes the propagation of the 
maximum up from the point Xq properly. 

For a propagation of a kink front that experiences a deformation, one can 
check the translation of the second derivative of the shape, keeping track of the 
turning-point of the kink shape (Fig.4). In this case the second order ("three- 
point" ) P V turns out to be the relevant velocity of propagation. 
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Figure 2: A damping deformation of a signal tp(x,t). A natural attribute allowing to 
describe propagation is the peak location (maximum) 



Summarized, the local concept is based on local attributes of the wave shape 
being traced in course of propagation. The phase of an arbitrary wave in some 
point is characterized by the local Taylor expansion in this point. The local 
attribute of N-th order is therefore the N-th spatial derivative of the wave shape. 
The phase velocity of N-th order (N-PV) is the velocity of propagation of the 
corresponding N-th order attribute. 

Finally, it should be noted, that the defmiton of phase velocities of order zero 
and one, the t>( ) and V(r> respectively, admits a straightforward generalization 
to two-, three-, and higher-dimensional propagation, while the phase velocities of 
second ( f(jn) and higher orders inherently contain a certain element of ambiguity 
in their definition, since a possibility to choose a second-order attribute to be 
traced is not unique is not unique [19]. 



4 Examples 

Let us consider the conventional 1+1-dimensional wave equation 



1 d 2 



d 2 



a 2 dt 2 dx 2 



ip(x,t) = 



possessing translational solutions of the form 

ij(x,t)=ij(t±-). 

(X 



(4.1) 



(4.2) 
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It is easy to check that in this case the PV's of all orders defined above are 
identical and read 

v (N) (x,t)=a, N = 0,1,2,... (4.3) 

which is nothing else but the classical phase velocity ±a, thereby satisfying 
the a priori definition of "phase velocity" itself as a medium constant in (4.1). 

Let a propagating shape \& now be subject to a temporal damping similar to 
(Fig.2) with 

t) = Mt - -)e- xt = V(0K A * (4.4) 
a 

that obeys the wave equation (a typical field equation with a dynamical dissipa- 
tion and a mass term) 

fid 2 d 2 2\ d A 2 | 

The ordinary 0-PV velocity reads 

"(0) = a (l " A|) , (4.6) 

where the prime denotes the derivative of ip with respect to its argument 
(j) = t — x / a. 
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Figure 4: The propagation of a growing kink ("tsunami model"). The turning point 
is chosen to trace the propagation 

This result provides a velocity with bad physical features: the velocity grows 
for a descending shape, for an ascending shape it decreases, can even be neg- 
ative, and it diverges exactly for the peak point (under the condition that the 
(measured) amplitude \1/ as well as parameter a, A have positive values). 

A relevant physical velocity in this case is for instance the 1-PV 

v {1) = a (l - X^j (4.7) 

providing for a peak being traced exactly the canonical phase velocity entering 
in the wave equation (4.5). 

For this shape further PV's of higher orders are given by (3.5): 

v (N) =a(l- = a (1 - AOogV"))- 1 ) (4.8) 

where ip^ denotes the n-th derivative of ip with respect to its argument as 
mentioned above. For an amplified signal as in Fig. 3 we can e.g. change the sign 
of A. Especially, for the case of a kink (Fig. 4) 

x 

ip(t ) = ip((f)) = arctan0, (4.9) 

a 

the spectrum of phase velocities reads by comparison : 



v (o) = a (l + A(l + 4> 2 ) arctan0) , 

V(i) = a (l - xl ~^f) » ( 4 - 10 ) 

The ordinary 0-PV has an oscillating sign at A and is multiple defined because 
of the arctan function . Therefore it cannot be interpreted as a well-defined 
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physical velocity. If the point being traced should be labelled by a derivative 
attribute, it turns out to be a physically inconvenient choice since the shape 
possesses no real peaks, whose propagation could be traced. Moreover the velocity 
t>(7) diverges at = 0. 

The possible labelled attribute is also the turning-point traced by the 2-PV. 
The velocity v^n) also possesses two singularities at = ±l/\/3 which do not 
coincide with the labelled point = 0, so it can be successfully followed up at 
the measurement. 



Historical remark 



The definition of a wave velocity proceeded with an extension of a translational 
solution of (4.1): 



if)(x,t) = if>(t - -) = if>[-(u;t- -x)],= Mut- kx), k = —, (4.11) 
a uj a a 

where the parameter cu has been interpreted as a frequency of a necessary perodic 
harmonic function 0. usually e ±tx as mentioned above in the introduction. It is 
not surprising that the 0-order PV provides in this case the value ou/k, canonically 
interpreted phase velocity of periodic wave 

In case of any interrelations between k and uj that are not encountered in 
the wave equation, in particular any dispersion relation between frequency and 
wave-number, the PV v (o) is in fact a proportionality factor 

duo = V(p)dk, (4-12) 

which is identical with the classical definiton of a group velocity [1]. 

At this point it should be recalled, that the original idea of the group velocity 
U, as summarized e.g. in [28], 

was of the similar form as the recent definition (3.2) of 0-PV. It should be noted 
here that equation (4.13) involves local derivatives of the non-local parameter, i.e. 
the wavelength A. For a variable wavelength much smaller than a vicinity of 
the point, this group velocity provides therefore a natural approximation to the 
zero-order phase velocity V(q). 

In the present approach this feature appears as a physical phase velocity 
following in a straightforward way from the interpretation of phase propagation 
and does not require an interpretation of uj and k as a frequency and wave number, 
as well as a constancy of some group respectively. 

It is noted in passing, that in the case of a kink there are no suitable definitions 
of a wave group and of a concerned group velocity for this solution since a Fourier 
decomposition does not work on a non-compact support. 
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5 Lorentz covariance 



The zero order PV (3.2 ) possesses an interesting property, i.e. a local covariance 
in the sense of special relativity, as shown now. This is not the case for PV's of 
higher orders. 

The PV f(o), measured in some stationary system X, takes in some other 
system X' , moving with a constant speed V, via the Lorentz transformations 

3 , = xpVt = i = ±IjL (5.1) 

the form 

v' (0) = M±^. (5.2) 

c 2 

This means that the zero order PV respects the relativistic velocity addition. 
Especially, a subluminal zero order PV remains also subluminal in any other 
moving system X' . 

This result should not be a surprise, since the definition of the i>( ) is implied 
by the condition: 

ip(x,t) = const, (5.3) 

which remains to be of the same form under arbitrary transormations 
x = x(x',t'); t = t(x',t f ), implying 

dip(x, t) = ^-dt + ^dx = (5.4) 
at ox 

for the first order differential form (or simply first differential), which possesses 
a form-invariance property under transformations. 
Since 

V{0) ( x,t ^ = ~dl ^'^ 

per definition, the equation (5.4) inplies the definition (3.2) of 0-PV, so it should 
behave under space-time transformations as a usual velocity of a matter point. 

The first order PV (3.4), by comparison, evaluated in some system X trans- 
formes in the moving system X' to: 



V (I) ( V 2 Q2 Q2 y Q2 -, , • (5-6) 

I c 4 dt 2 ~T dx 2 Z c 2 dtdx J V 



If the pulse if)(x,t) obeys the free wave equation (4.1), the transformation 
becomes 
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(l + £ ) „<„ + 2V 



= " (1+5)+^ ' (5 - 7) 

i.e. a relation that should be called the "first order velocity addition". It differs 
obviously from the corresponding transformation of t>( ), since the definition of 
the velocity results from the condition 

whose form is explicitely non-invariant under space-time transformations. 

It can be shown that a subluminal first order PV is still preserved by this 
transformation as well. Note, that for a signal which does not obey the free 
equation (4.1), this restriction is in general not guaranteed anymore. 



6 A global velocity of signal transmission and 
"dynamic separation" 

The discussion of local velocities was aimed towards an evaluation of global fea- 
tures of signal propagation, namely the propagation through a finite spatial inter- 
val during a finite temporal interval. In other words, a global velocity, practically 
measured, means roughly the length of the distance Ax traveled by a traced 
attribute divided by the time interval At. 

The local PV's of N-th order analyzed above can be interpreted as a first 
order differential equations of the form 

dx 

V(N)(x,t) = —, (6.1) 

that can be illustrated graphically as a field of isoclines (Fig. 5). Here the local 
PV is the tangent function of the tangent vector of the isocline, and the averaged 
(total) velocity between (to,xo) and (ti,xi) is represented by the tangent of the 
hypotenuse of the triangle {(to,xo), (ti,xi), (to,xi)}, (cf. Fig. 5). 

This procedure is elucidated best by some clear and well known examples. 
Consider the propagation of a translation mode of the form 

^(x,t)=^(#-MaO):=W) (6.2) 

describing, for instance, an electromagnetic field propagation in an inhomoge- 
neous dielectric medium. In this case, if the transversality of the wave not as- 
sured, the differential equation for any components of the field can possess e.g. a 
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Figure 5: Phase velocities as a family of isoclines v^(x,t) and averaged global veloc- 
ities between two events (measurements) 

form like this: 

• h.V.r)^r <> (6-3) 



c 



.2 



dt 2 dx 2 dx 



where n(x) is the spatially dependent index of refraction. From the Ansatz (6.2) 
it results for the k{x) from the optical inhomogeneity over the spatial coordinate 

£ C J 

Then the local zero order PV is provided by eq.(3.2) and results in the first order 
equation 

v ^ t) ^ = WY (6 ' 5) 

that is a quite trivial result. 

We proceed now with phase velocities of higher orders. For the 1-PV velocity 
t>(7) defined by eq.(3.4) the same procedure leads to 

at n[x) - ||t7 In n[x) 

where x and t denote the space and time location of the first oder attribute being 
traced; ip'andip" are the first and the second derivative of the wave function if) 
with respect to its argument (cf.6.2). Similarly, for the 2-PV we obtain: 

V(ii){x,t) = — = —— r , C 2 n „ j? ■ (6-7) 

dt n{x)-3 I ^ T7 lnn{x) + ^^f TT 
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Suppose, we measured at the moment to some value of the fisrt order attribute 
ip'(x ) in the point x , and we are interested, where is this attribute to recover 
at the next moment t\. For a maximum or a minimum point the result does not 
distinguish from the zero order velocity, as in the case considered above. Now, 
let we measured at t some values: 

V>'(0, x ) := ij' and ^"(0, x ) ■= <; A := V'(0, x )/< (6.8) 
An explicit integration of eq.(6.6) leads to 

cfa - to) = J n(x)dx - A| In ^ (6.9) 

for the interval of two events, where the signal attribute choosen for tracing 
has been checked at the time to in the point x and afterwards in x\ at t±. Let us 
assume that the point (to,xo) was not a critical point (i.e. of the knot type) of 
the equation (6.1). Then the global transition velocity between two points xq,xi 
is evaluated unique in such a way: divided by the spatial interval x\ — x the 
eq.(6.9 ) provides the factor 7: 



ti — to 

7 := c- 



X\ — Xo X\— Xo 



X\ 



J n(x)dx - A| In 



n(x ) 



LEO 



(6.10) 



It follows, that a value of 7 less than 1 is not prohibited in principle even for 
real n > 1, that means a superluminal averaged global velocity along the finite 
distance. A similar result appears for 2-PV and higher orders. A numerical 
examples for dielectric media of special modulation confirms this possibility [27]. 
The discussion of physical interpretation of this result is still open. 

It is also worth to notice, that for the first order PV (as well as the second 
and higher orders) in media with a variable refraction index n(x) an essential 
dependence on the parameter £ enters. 

The meaning of £ can be derived from the given form of the solution ip. It 
can be interpreted, for instance, as a frequency factor for a periodic mode or 
a damping factor for an evanescent mode. This parameter simply means an 
enumeration of solutions of a solution family (space of solutions). The merit 
gained above is the separation of these solutions with respect to the parameter £ 
by the different first (and higher)-order PV's through the inhomogeneity of the 
medium. 

In the case of a periodic wave for instance, it can be interpreted as a dispersion 
although an explicit frequency- dependence of n(x) was not assumed. Moreover, 
the parameter £ has not necessarily to be interpreted as a frequency of some 
temporally periodic oscillation, but rather as a time component of the space- 
time wave vector {u, k} for a translation mode (4.2) in a general form. The 



15 



phenomenon considered thus has another origin and is much more general as a 
conventional non-localized frequency dispersion n(uj) 

Thus we established for the first (and higher) order PV's in optical inhomoge- 
neous media the essential dependence on the t-component u of the wave transla- 
tion vector, especially frequency, even for local non dispersive media, which could 
be called " dynamic separation" . The global averaged dynamic separation of eq. 
(6.9) survives as a corollary of the local separation of (6.6). This phenomena does 
not occur for the ordinary zero order PV. 

7 Concluding remarks 

The inherent inconsistency of the classical subjects of phase velocity and group 
velocity, (as well as signal velocity based therein) has been discussed. The inap- 
plicability of these concepts for actual studies in photonics, near-field and nano- 
optics has been shown to result from the essential non-locality of these definitions. 

An alternative approach for description of a propagation velocity has been 
proposed. It is strictly local and is based on the natural assumption of an ordinary 
measurement procedure. It does not need any a priori condition of periodicity, 
frequency, groups and packets or other canonical attributes. 

The definitions presented above are applicable in a natural way for arbitrary 
waves and pulses. In a mathematical sense they describe a propagation of a 
perturbation in any field in space-time. Examples could be: an acoustic wave as 
a pressure perturbation, a gravitation wave on a fluid surface, a spin wave in a 
solid state a.s.f. In particular, the formalism is very suited for the description 
of particle propagation in field theory, where particles are considered as field 
perturbations. 

This approach results in the set of measurable propagation velocities. In 
zeroth order the propagation velocity coincides with the ordinary phase velocity 
and appears to be ordinary Lorentz covariant; further application gives rise to 
generate " Lorentz covariance of higher orders" . 

For propagation in inhomogeneous media, for instance for light in a medium 
with a space-dependent index of refraction, phase velocities of higher orders show 
some remarkable properties: 

a) they are not restricted in general by the speed of light in vacuum; 

b) they exhibit an essential dependence on the time component of the wave 
vector solely as a result of inhomogeneity, treated canonically as a " dispersion" . 
This appears to be a more general phenomena, and does not presuppose any 
periodic frequency and dispersive properties of media. It should be called for this 
reason "dynamic separation" [27]. 

Further applications, like photonics and optoelectronics are dealing with sig- 
nals i.e. non-localized objects, represented as a cut of wave shape between la- 
belled points possessing a certain spatial length (or temporal duration). These 
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two points are marked by a value of zeroth order attribute (e.g. a voltage am- 
plitude), or by a first order one (a block between an ascent and a descent of 
amplitude). In this sense the signal is considered as not a local, but at least a 
bi-local object. Its propagation can be traced by an averaged bilocal velocity, 
constructed from local ones in a natural way. 

The same concerns a propagation of single photons, especially of low energy, 
which are non-localized objects. For the case of evanescent modes we are dealing 
yet with a sufficient derealization. 

Another approach to the signal transmission is based on a pure local signals, 
which can be, for instance, a point of discontinuity, propagating away from the 
perturbation [18]. For these problems the velocities of higher orders, especially 
1-PV and 2-PV provide a suitable tool, as shown above. In this context are the 
velocities 1-PV and 2-PV of special interest, being treated in pure physical sense 
as a true speed of interaction. 
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